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Abstract 

It is shown that the low energy coefficients of the next-to-leading order (NLO) chiral perturbation 
theory needed to determine AI = 1/2, K vrvr decay amplitudes can be fixed by calculating 
Ktt TT amplitudes on lattice. Unlike using NLO K vrvr amplitudes proposed by Laiho and Soni, 
simulating Kir vr transitions on lattice does not require evaluations of s-channel disconnected 
diagrams which have been an obstacle in practice. 
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I. INTRODUCTION 



Chiral perturbation theory(xPT) has been a useful tool in understanding the physics 
of mesons. In particular, it has been employed to extrapolate important quantities to the 
physical pion mass from lattice calculations performed using somewhat heavier mesons. 

Employing leading order chiral perturbation theory in the lattice QCD calculations of 
nonleptonic kaon decays was proposed by Bernard et al in Ref.j^. In this proposal, the 
effective weak operators are rewritten in terms of meson fields. The coefficients of those 
operators(low energy coefficients) are determined through lattice simulations where lattice 
calculations of ^ vac and K ^ ir can be used to determine the more difficult K ttvt 
amplitudes. This calculation has been done [sl. \l4 1 with quenched ensembles. However, the 
results did not show a good agreement with experimental observation. 

Since then computing power has drastically improved and we can afford to generate 
unquenched lattice ensembles. The chiral perturbation theory treatment has also been 
extended with the next-to- leading order calculation done by Laiho and Soni|l^. However, 
at next leading order K vr7r/=o calculations in addition to K ^ n and K vac on lattice 
are needed to determine all K tttt amplitudes with physical kinematics. 

Unfortunately, there are significant difficulties in lattice calculations of K 7r7r/=o tran- 
sitions. Those difficulties basically result from the existence of s-channel disconnected di- 
agrams. The numerical evaluation of those diagrams turn out to be very hard because of 
the exponential decay of signal to noise ratio. In order to avoid the difficulties, the author 
have proposed to use Kn ir transitions in Ref. 9|. The most notable advantage of using 
Kit 71 amplitudes is the elimination of s-channel disconnected diagrams. Since K^t[~ 
state is the highest weight state of / = 3/2 states, the Iz cannot be zero after the interaction 
with weak operators whose = —1/2. This non-zero Iz guarantees the absence of s-channel 
disconnected diagrams. A more detailed discussion about the difficulties can be found in 
Reffl. 

In this paper, the detailed xPT formulas for Ktt — > vr processes are presented. In the next 
two sections, we fix notation by expressing the effective weak operators in QCD as elements 
of definite irreducible representation of the SU{?>)l ® SU{?>)r chiral symmetry group and 
their corresponding expressions in terms of meson fields follow. Then, it will be explicitly 
shown that the LECs needed to determine the physical K vrvr decay amplitude at NLO in 
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xPT can be determined by analyzing the chiral expansion of Ktt — ^ vr amplitudes at simple 
kinematic points. Those effects of finite volume are discussed which could be an issue when 
two particle states with non-zero relative momentum are considered. 

II. EFFECTIVE WEAK OPERATORS 

The operator product expansion(OPE) can be used to express K — >■ tttt decay amplitudes 
in terms of matrix elements of the AS* = 1 weak effective Hamiltonian, 

{7rn\nAs=i\K) = ^V^,V:,J2^^if^){^^mf^m. (1) 

In essence, the OPE separates two important physical scales : the Cj(/i)s, called Wilson 
coefficients, which contain the short distance physics which can be calculated by QCD 
and electroweak perturbative techniques and the matrix elements, {7ni\Qi{fi)\K) which are 
determined by the long distance physics for which nonperturbative methods are required. 
A thorough discussion can be found e.g. in Ref. ^. 

The Qi{fJ.) are four quark operators of which there are 10, 

Qi = Sa^^.{l - ^^)uaUbl''il - 7^)db, (2) 

Q2 = Sa7M(l - I'^Ubl^il - l')da, (3) 

(4) 
(5) 
(6) 
(7) 

(8) 
(9) 
(10) 
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Qio = 2^a7M(l - 7^)4 MftT^'ll - 7^)ga. (11) 

The Qi and Q2 are called current-current weak operators. The operators Q3 — Qq arise from 
QCD penguin diagrams whereas Qj — Qio are from electroweak penguin diagrams. 

In order to use chiral perturbation theory, these operators must be written in terms of 
elements of irreducible representations of the chiral symmetry group. The relevant four- 
quark operators can be arranged into irreducible representations of the SU{3)l ® SU{3)r 
chiral group with definite isospin as follows : 



y{3/2) 

^27,1 ~ 


{sd)L [iuu)L - idd)L] + {su)L{ud)L, 


(12) 


v(l/2) _ 
^^27,1 ~ 


{sd)L [{uu)l + 2(J(i)i - 3(ss)l] + {su)L{ud)L, 


(13) 


v{l/2) _ 
^8,1 - 


{sd)L{uu)L - {su)L{ud)L, 


(14) 


i>{l/2) 
^8,1 - 


{sd)L [{uu)l + 2{dd)L + 2{ss)l\ + {su)L{ud)L, 


(15) 


v(V2) _ 

-^8,1 ~ 


{sd)L [{uu)r + {dd)R + {ss)r\ , 


(16) 


V(V2)C _ 
-^8,1 — 


{ {sd)L [{uu)r + {dd)R + {ss)r\ Y , 


(17) 


^;(3/2) 
->^8,8 — 


{sd)L [{uu)r - {dd)R\ + {su)L{ud)R, 


(18) 


^,(3/2)c _ 
-•^8,8 — 


{{sd)L [{uu)r - {dd)R] Y + {{su)L{ud)RY , 


(19) 


V(V2) 
->^8,8 — 


{sd)L [{uu)r - {ss)r] - {su)L{ud)R, 


(20) 


V(V2)C _ 
-^8,8 ~ 


{{sd)L [iuu)R - iss)R]y - {{su)Liud)RY . 


(21) 



where we follow the notation of Ref.[lJ]. Operators are classified by the Lorentz structure 
L ® L and L ® R and represented by X and y respectively. The irreducible representation 
to which the operator belongs is given in the subscript and the short hand notation 

{sd)L = r^,{l - l5)d'' , {sd)R = s\{l + -f,)d' (22) 

is used. Because of invariance under the Fierz transformation, there are no color mixed L^L 
operator appearing, but for L ® R operators, their color mixed versions are represented by 
superscript c : 

{{sd)L{sd)RY = s%(l - 75)^?" Am(1 + 75)c^^ (23) 
Finally, the isospin of the operators is also given in the superscript. 
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In terms of this basis, the four-quark operators are rewritten as 

^ _ 1 y(\l2) 1 y(\l2) 1 ;p(l/2) 1 ;p(3/2) .^ 

VI - 2 8,1 + Yo^'^ 15 3^27,1 5 

^ _ 1 ;^(l/2) 1 ;p(l/2) 1 ;p(l/2) 1 ;p(3/2) 

V2 - ~2 10 15 3^27,1 5 l^^j 

Qs = ^4Y^^ + ^^ir> (26) 

Q4 = Q2 + g3-Qi, (27) 

Q5 = 3^^;f\ (28) 

Q6 = 3^^;f ^ ^ (29) 

Q7 = ^[3^Sf +3^ir]> (30) 

Qs = ^[3^^f ^ + (31) 

= \Qi - \Qz. (32) 

Qio = Q2-\Qz^\Qx. (33) 



The strategy is to measure the matrix elements of the X and y operators with quarks 
somewhat heavier than physical quarks and to use the chiral expansion to extrapolate the 
matrix elements to the physical quark mass. Then, we can recover the matrix elements of 
Qi(/x)s from Eqs.dMHSSD- 



III. CHIRAL PERTURBATION THEORY 



Chiral perturbation theory (xPT) is based on an effective field theory of the low energy 
sector of QCD. Its fundamental degrees of freedom are the lowest mass pseudoscalar mesons 
which are the Goldstone-bosons arising from spontaneous chiral symmetry breaking. Be- 
cause of the non-linear transformation of the Goldstone-bosons under the symmetry group, 
the meson fields appear in the field S given by: 



exp 



/ 



(34) 



where S belongs to the (3,3) representation of the SU{3)l ® SU{3)r, the 3x3 matrices 
t°- are proportional to the Gell-Mann matrices with tr{tatb) = Sab, and the 0" are the real 
pseudoscalar-meson fields. The quantity / is the meson decay constant in the chiral limit, 
with equal to 130 MeV in this notation |lo| . 
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The leading order (0(p^)) strong Lagrangian is given by 

4? = ytr[9,St9^S] + ^trix^S + St^], (35) 

where x is a diagonal mass matrix with its diagonal elements (m^, m^, m^). They are related 
to the meson mass via 



^1,+ m^o , X 

Bo = = = (36) 

rriu + md mu + rUs ma + 



The leading order weak Lagrangian is given by 

= aitr[A69^St9^S] + a225otr[A6(x^S + Stx)] 

+a274{^^d^j:)'^{j:^d^^ + asstrlAeStgS] + H.c. (37) 

where terms with coefficients ai and a2 belong to the (8, 1) representation of the SU{3)l ® 
SU{3)ji group and the last two terms belong to the (27,1) and the (8,8) representation 
respectively. 

The next-to-leading order(0(p^)) weak operators contributing to kaon decays are given 
in Ref. 6|, with the effective Lagrangian, 

£<j) = j2 '^'or + E + E =<o!''". (38) 



The (8, 8) and (27, 1) operators will not be considered in this paper since t 
operators can be determined without AI = 1/2, K tttt simulations 
forms of the (8, 1) operators are : 



le LECs for those 
The explicit 



(^(8,1) 




(^(8,1) 


= tr[A65]tr[5], 




= tr[X,P% 


0(8,1) 


= tr[A6P]tr[P], 




= tT[Xe[S,P]], 




= tr[X,{S,L'}], 




= tT[XeL,,SL^], 


A^(8,l) 


= tr[A6L,,]tr[{L^5}] 




= tr[A65]tr[L2], 




= tr[A6[P,L2]], 




= tT[x,{L,,d,wn], 




= tilXeW^^W'^'^] 



with S = 2Po(x^S + Stx), P = 2Po(x^S - Stx), = lE^d^E , 1^^^ = 2{d^L, + d,L^), 
and {Xe)ij = 63i62j. 

Finally, the next-to-leading order strong Lagrangian relevant for kaon decay amplitudes 

is 

4? = E^^^1"^ (39) 
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and the explicit forms for the are 





= tr[L2]2, 


of) 


= tr[L^L,]tr[L'^L- 


of) 


= tr[L'L% 


of) 


= ti[L^]ti[S], 


of) 


= tT[L^S], 




= tr[5]^ 


^8 









(40) 



IV. ROLE OF Ktt^tt MATRIX ELEMENTS 



In this section, we will show explicitly that those LECs necessary to calculate the physical 
K ^ iTTT amplitude can be determined from K vac, K n and Ktt vr amplitudes. 

There are five types of operators depending on the representation of the chiral group 
and isospin : (27,1) AJ = 3/2, (27,1) AJ = 1/2,(8,8) AJ = 1/2,(8,8) AJ = 3/2 and (8,1) 
A/ = 1/2. In this paper, only (8,1) A/ = 1/2 operators are discussed. LECs for other 
operators (27,l)/(8,8) A/ = 1/2 or 3/2 can be determined without the use of K 11111=0 
amplitudes |l5[|. 

In order to check that all the LECs sufficient to reconstruct the physical K irir ampli- 
tude can be determined, it is enough to look at the analytic terms. Since matrix elements, 
meson masses and the meson momenta are calculated on lattice, if we insert those quan- 
tities into the following xPT formulae a set of linear equations of LECs will be obtained. 
Presumably, the contribution of the logarithm terms would not render the linear equations 
singular. So, we omit the logarithmic terms in the following formulae. 

Although K vac and K ir calculations have already been done, for example, in 
Ref . [lo| . we present the formulae for completeness. The analytic terms in the K — > vac 
amplitude come from diagrams Tl and T2 in Figil] and they are 

/Q|^(8.1) I j^O^ - ^'^^ ^"^^ ~ ^"^^ ~ "^^^ ^ ^ ^ ^^"^"^ (41) 

From this calculation, one can determine a2, 62, ei + 65. If CPS symmetry|l] is realized, 
rriK = ''tItt; the K vac matrix elements vanish. Thus, non-degenerate quark masses must 
be used for the determination of these LECs. 

The diagrams El and E2 in Fig I2] generate the analytic terms in the K ^ it amplitudes, 
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which are 

+ ^2 ( - 16 (ei + 62 + 65) 

+ (16 (eio - 635) {pk ■ K) -8 (62 + 263 - 265) ml) 



m 



K 



+64639 (px -Kf 

+8 (611 -2635) (pa- -M^'). (42) 



By varying momenta and masses, one can determine ai, 639, 611 — 2635, 610 — 635 and 63 — 65 
when combined with K vac results. 

For A'tt 71, there are five diagrams: Al, A2 and CI in Fig|T]and Gl and G2 in FigS 
The analytic terms are 



g 

-a2 {m\ - ml) - Aai{{pK ■ K) - {Pn ■ K)) 

32 X6 
+ Y (ei + 62 + 65) m\^ + Y + + 363 - 265) mlm\ 

-16 ((610 - 635) {pk ■ K) + (2613 + eis) (Ptt ■ K)) m^K 
16 

— (36i + 762 + 363) ml - 64639((Pa' • - {Pn " Kf) 

-8 (611 - 2 (615 + 635)) {pk ■ K)ml 

+8 (2610 + en + 46i3 - 4635) {p^ ■ K)ml 

+ 16635(px ■P-k){{PK ■ K) - (Ptt ■ K))), (43) 



(7r-(A;.)|0(«'i)|ir°(pA.)vr-(P.)>^n.^2 = 

4i {-ml + 2{pK ■ Pn) + {pk ■ K) + (Ptt • K)) {ml - m|-) (eta + 4 (61 + 62 + 65) m|- + 262m^) 

3/3 (m2 + (pj^ ■ p^) - {pk ■ K) - {p^ ■ K)) 

(44) 
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and 

il6a2{m1^ — ml) 



-12k^ ■ pkK ■ Pn (2/^1 + L2 + L3) 

+6pK ■ Pn {K ■ Pn (4Li + 2L2 + Ls) + K-PK (4L2 + Lg)) 

+2K ■ Pn (12Li + 3L3 - 8L4 - L5) - ApK ■ p^ (2L4 + L5) 

- (2L4 + L5 - 2 (2L6 + Lg)) - 2K ■ Pk (2L4 + L5) m^, 
+2A;^ ■ Pk (6L2 + 3L3 + 5L4) - 2A;^ ■ p^L4m^ 

- (IIL4 + 4L5 - 6 (SLe + 2^8)) mlml - 2pK ■ p^ (6L2 + 8L4 + 3L5) ml 
+ (L4-L5 + 2(L6 + L8))mf 



where the Li are Gasser-Leutwyler coefficients and it is assumed that those coefficients are 
aheady known. 

Since the analytic terms of Kit n amphtudes are quite complex, we isolate the new 
coefficients which must be determined from Kn tt amplitudes by inserting LECs that can 
be computed from K — > vac and K ^ tt amplitudes. One can easily see that only three 
coefficients remain to be determined which are from the G2 contribution, 



I 

16635 (Pi^ ■ Pi^){{PK ■ K) - (Ptt ■ K)) 



-32ei3(p^ ■ K) [m]^ - ml) 

-16ei5 ((Ptt ■ - {pk ■ K)ml) j . (45) 

In order to show explicitly that the needed coefficients can be determined, we choose kine- 
matics where the initial kaon and pion are nearly at rest while the final pion has momentum 
27r/L. With these kinematic points, 

{^^{K)\0^'^'^K\pK)ir~{p^))^^^^^^^ = 

mW ((635 - e^^){W^ml - WKml) - e^^{Ew + W^){ml - ml)) 

P 

(46) 

where W is the energy of final pion and Wk and are the energies of initial kaon and pion 
respectively. Although very small, we take into account the momentum of the initial parti- 
cles. When Kti tt transitions are simulated on lattice, this small momentum originates 
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from interactions between the kaon and the pion and we have no control over the direction 
of this momentum. An average over sohd angle must be taken. More discussion of this point 
is given in the section |Vl This effect is already included in the above formula. Now, one can 
see that 635 — 615 and 613 can be determined by varying mx and rrij^. Then, combined with 
K — ^ vac and K ^ tt results, one can determine the physical K irir matrix element : 

/ 2 2\ 

— (m^ - m^) ■ 

(^(eio - 2ei3 - ei5)m|- 

+ (-2ei - 462 - 263 + 2610 + en + 46i3 - 4635 + 8639) j . 

Note that the SU(3) limit cannot be used because Eq.( l46l) vanishes in that limit. More- 
over, there is a further restriction on the choice of momenta if one considers more general 
kinematics. The tadpole contributionfEq. (1441) ) can diverge for some choices of momenta. 
This is a disadvantage of using transitions with an unavoidable energy- momentum injection 
such as the current proposal. 



V. FINITE VOLUME EFFECTS 



In the previous section, it is shown that the unknown LECs can be determined by mea- 
suring weak matrix elements from lattice calculations with mesons nearly at rest. We choose 
to use a minimal set of matrix elements for the sake of proving sufficiency of the calculation 
of only these types of amplitudes. However, it would be useful to measure more matrix ele- 
ments rather than the minimal set so that one can reduce statistical errors in determining 
the LECs and test for the consistency of the xPT expansions. 

An obvious approach is to include matrix elements with mesons of different masses, but 
this requires new sets of configuration if we want to use unquenched chiral perturbation 
theory. A method which does not require new ensembles is to use mesons with non-zero 
relative momenta. In this case, the finite volume effects on the resulting matrix elements 
are not exponentially small. This has been studied by Liischer and Lellouch in Ref. 11 1 
and it is generalized to the case where the total momentum is not zero in Refs.j^, Q]. This 
generalization is quite useful since using a system with non-zero total momentum is the only 
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method, at the moment, of creating two particle states with non-zero relative momentum 
without generating new gauge ensembles. 

One can also consider matrix elements whose final meson states have non-zero momentum. 
In such cases, the result of Refs.js, 8] has to be generalized further. This is straightforward 
and the derivation is given in appendix [Bl The final result is 



(7r(k)|i7^(0)|ir7r,(E,P)) 
1 1 



TT 



d^E* (7r(k*)|i7H^(0)|i^(q)7r(-q)) 



(47) 



where the left-hand side represents a matrix element in a finite box while the integrand on 
the right-hand side represents one in an infinite volume. Here, E is the energy of the K% 
state which can be measured from lattice calculations. Similarly, P is the total momentum 
of the Kti state which is imposed explicitly by the operator creating the state. Finally, k is 
the momentum of the final pion. The starred variables are those Lorentz-transformed into 
the CM frame. In particular, 

E* = VE^ - P2 (48) 
which in turn gives the Lorentz transformation angle j3 = P/E and the q* is defined from 



E* 



*2 



(49) 



The k* is the result of the Lorentz transformation with the above (3 of the four-momentum 



of the final pion state, {^/to^ + k^, k). py is a function of g* and its definition is given in 

states. 



Refs. js], Q] . Roughly, it can be interpreted as the density;^o 

In Eg. (1471) . the difference from the result of Refs. 
the integration over solid angle. Since the final meson state has non-zero momentum and 
cannot be used to constrain the initial state to be an S-wave, this must be done explicitly. As 



is the explicit appearance of 



such, the initial state in the matrix element (7r(k)|ifvy(0)|-R'vr, {E, P)) must be in an S-wave 
state [l^, as seen from the CM frame. This means that the S-wave state must be explicitly 
generated in the lattice simulation. Fortunately, S-wave states are the lowest energy states 
which can appear, so we can generate the S-wave Kn state by using any operator which 
has an overlap with it and evaluating its correlation function at a large Euclidean time 
separation limit. 
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D2 D3 D4 



FIG. 1: Diagrams involved with tadpoles. T1,T2 and T3 are relevant to K ^ vac amplitudes while 
others are for Ktt vr amplitudes. Gray(black) blobs are weak(strong) vertices. Large blobs are 
the leading order terms while the small blobs are the next-to-leading order ones. 

VI. CONCLUSION 

In this paper, it is shown that the LECs necessary for the NLO xPT calculation of 
the physical (^mi\0^^'^^K'^ matrix elements can be obtained from lattice calculations of 
Kit — > 7r,K ir and K — * vac processes. The important point is that simulations of 
K — > TTTTi^o transitions, which are very difficult, can be avoided. 

Although we establish this result by using a minimal set of kinematics, the finite vol- 
ume effects on K ^ -Kir matrix elements are also discussed, effects which should be taken 
into account if one wants to improve the statistics and the control of systematic errors by 
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o 

El 



-<y- 
E2 



Fl 




F2 






HI 




H2 




H3 



FIG. 2: Diagrams not involved with tadpoles. E1,E2 and Fl are relevant to K ^ n amplitudes 
while others are for Kir — > vr amplitudes. Gray(black) blobs are weak(strong) vertices. Large blobs 
are the leading order terms while the small blobs are the next-to-leading order ones. 

including more kinematic points. In particular, the case where the final one particle state 
has non-zero momentum is discussed, which requires a slight generalization of the results of 
Ref.y. 

As discussed in Ref.jsl, using Kir tx processes allows us to avoid technical difficulties 
such as s-channel disconnected diagrams. However, the mixing with lower dimensional 
operators is still present. One may apply the subtraction scheme used in Ref.jsl. 

The current computing resources are powerful enough to simulate Kn tx transitions. 
With these NLO xPT formulae, we can improve the calculation of — ^ vrvr matrix elements, 
which in turn will allow us to evaluate e'/ e more accurately. 
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APPENDIX A: LOGARITHMIC TERMS 

In this section, the logarithmic terms (diagram A3,F1,F2,D1,D2,D3,H1,H2,H3) are given. 
1. The Integrals 

In order to write the transition amplitude in a simpler form, the following notations are 
used : 
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1 




1 




{2t:Y {q 




— mf (g + ky 


- ml 








1 




{27ir {q 


-kf 


— m\{q + ky 


- ml 


d'^q 






qv 




{27rY {q 


-kf 


- mf (g H 


vky 


-ml 


d'^q 


q' 




1 




{2tiY {q 


-ky 


- (g H 


vky 


-ml 


d'^q 


qi^ 




q' 




{2tiY {q 


-kf 


- (g H 


vky 


-ml 


d'^q 


q' 




q' 




{2nY {q 


-kf 


- (g H 


vky 


-ml 



Bo{mi,m2,k) = /i' 
Bof,{mi,m2,k) = fx' 
Bo^^{mi,m2,k) = fi" 

Bi{mi,m2, k) = /i' 
Bif,{mi,m2,k) = /i' 

B2{mi,m2, k) = fi" 

Furthermore, when there is a Lorentz contraction, we use the following abbreviation 

k-Bo = Bo^k^^ 
k-B^ = B.^k^' 
p-Bo-k= p^Bo^uk" 
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2. Diagram Hla 

Here, g = - K), 

AaiK ■ Bo (g, mx, mx) ■ Pk + SaiK ■ Bq (g, m^, m^) ■ px + 4aip„ ■ Bq (g, uik, mx) ■ Pk 

2 

+8aip^ ■ Bo (g, m^, m^) ■ px - -«2 (4m|- - 4m^) 5i (g, mj^, tuk) 

+ (-2m|^ - 5m^ + ■ Pn) {-ml - K ■ px + K ■ Pn + Pk ■ Pn) Bo (g, mx, mx) 

+ ^^2 (-2m|. - 5m^ + ■ Pn) (4m^ - 4m^) 5o (g, m^, mx) 

+ ^^2 (8m|. - 8m^) (llm^ - lO/c^ ■ p^) Bq (g, m^, m^) 

+ (-6m^ - 6A;^ ■ px + 6A;^ ■ + 6px ■ Pn) (H"^^ - 10A;,r ■ I^tt) (g, m^, m„) 

2 4 
+ -^2 (4m|- - 4m^) Bo (g, m^, m,,) + -02 (8m^ - 8m^) Bi (g, m^, m^) 

2 

+ -aiml {3ml + • px - 3K ■ Pn - ^Pk ■ Pn) Bo (g, m^, m^) 
y 
2 

(2m| + 6m^ + /c,r ■ Pi^ - GK ■ Pn - Pk ■ Pn) Bi (g, mx, mx) 

4 

+ -«! (Sym^ + 6k^ ■ Px - 36k^ ■ p^ - 6px ■ Pn) (g, m^, m^) 

4 4 16 \ 

--aiBi (g, m,^, m^) ml + -aiB2 (g, mj^, mi^-) - —aiB2 (g, m^, m^) ) 
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3. Diagram Hlb 

Here, q = ^{p^ - k^), 

-SaiK ■ PkK ■ Bq (g, m^, m^) + 3aiK ■ PnK ■ Bq (g, itik, m^) - aiK ■ PkK ■ Bq (g, m^, m^) 
+aiK ■ PnK ■ Bo (g, m^, uik) - QaiK ■ Bi (g, tjik, m^) - 2aiK ■ Bi (g, m^, m^) 
-SaiK ■ Bo (g, m^^, m^) pk ■ Pn - aiK ■ Bo (g, m^, m^) Pk ■ Pn - SaiK ■ PkPk ■ Bo (g, m^, m^r) 
+3ai/c^ ■ p^p/^ ■ Bo (g, m^^, m^) - "iaiPK ■ P^Pk ■ Bo (g, trk, m^) - aiK ■ PkPk ■ Bo (g, m^, tjik) 

+aiK ■ PttPk ■ Bo (g, m^, tuk) - a^PK ■ P^Pk ■ Bo (g, m^, mj^) - Gaipx ■ Bi (g, m^, m^) 

5 25 
-2aipK ■ Bi (g, m^, m^i-) - -aiK ■ PkP-k ■ Bo (g, niK, m^) + —aiK ■ PkP-k ■ Bo (g, m^, m^) 

3 2 10 

+-aiK ■ Bo (g, m^-, m^) m|. - -aip^ ■ Bi (g, m^-, m^) + —aip^ ■ Bi (g, m^, mx) 

8 13 
+ -a2A;7r ■ 5o (g, m^, m^) m|. + -aiK ■ Bq (g, m^, rriK) m\ + 2*^1^''^^ ' ^0 '^'^) 



1 ^ . ^ o 5 



+ ^a2PK ■ Bo (g, mii-, m^) m^^ + -aiPi^- ■ Bo (g, m^, m^) - -aip^ ■ Bo (g, niK, 

-T^^iPtt ■ 5o (g, m^, m/<) + \aiK ■ Bo (g, m^, m^) ml - ^a2K ■ Bo (g, rriK, m^) ml 

13 8 
+ -aiK ■ Bo (g, m^, mi^) + -aipK ■ Bo (g, mj^, m^) ml - -a2PK ■ Bo (g, mx, m^) ml 

15 5 
+ o"iP/^ ■ ^0 (g, rur,, mx) ml - -aip^ ■ Bo (g, mx, m^) ml - —aiPn ■ Bo (g, m^, mx) ml 
z o io 

-2aiK ■ Bo (g, m^^, m^) ■ - lOaiK ■ Bq (g, m^, mi^) ■ p^ - 2aip^ ■ Bo (g, mx, m^) ■ px 

+^"1 {-m\ -ml + 2K ■ Pk - 2K ■ Pn + 2px ■ Pn) (SOA;^ ■ px - 23 (mj^ + ml)) Bo (g, mx, m^) 
2 

+ n«2 {mx - ml) (7 (m^ + ml) - 10/;;^ ■ pi^) Bo (g, mj^, m^) 

+ - 2A;^ ■ px + 2K ■ Pn - 2px ■ Pn) + ml - 30K ■ Pk) Bo (g, m^, mx) 

+^ai [-lAm\ - lAml + 30A;^ ■ px - ■ Pn + ^Pk ■ Pn) Bi (g, mx, m^) 

(-10m|- - 2ml + IS^^r ■ Pk - ■ Pn + ^Pk ■ Pn) {q, m^, mx) 

8 10 2 

-n"2 ('Tii' - "2^) Bi (g, mx, m^) + — ai^a (g, m^^, m^) + -tti^a (g, m^, m^^) 
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4. Diagram H2 

Here, q = -Kpk+Ptt), 

4 16 40 

-aiK ■ Bo (g, TTiK, m^) m\ - —aiK ■ Bi (g, m^, m^) + y ai/c^ ■ Bq (g, m^, m^) ■ Ptt 
+ ^aiK ■ Bo (g, m^:, m^) - (m|- - ml) {m\ + + lOpi^ ■ p^) 5o (g, m^, m^) 
+ (m^ - ml) Bi (g, m/^, m^) j 

5. Diagram H3 



//3 /^TT^ 
64 

— (-35aiA;^ ■ px + ^aiPK ■ Pn + 2 {25aik^ ■ p^ + 7a2 {m\ - ml))) Aq {mx) 
32 

— — [-Aa2m\ + Aa2ml + lOai/c^ ■ px - Tai/c^ ■ + SaiPi^ ■ p^) Ao (m^) 
32 

+— {^a2m\ - 8a2ml - 90ai/c^ ■ px + 4:5aik^ ■ p^, + 45aiPi^ ■ p^r) ("^r?) 

224 32 \ 

-64aiyli (mx) + — ai^i (m^) - y ai^i (m^) j 



6. Diagram D12a 

Here, g = Upx + Ptt) 



oi2{m\ - '^D 



(m^ - 3m^ + 4/c^ ■ + 4/c,r ■ Ptt + QpK " Ptt) + + ^^^Pk " Ptt) (g, mx, m„) 

2 

+ n ("^A' -^1 + '^f^TT ■ PK + '^K ■ Pn + 8px ■ Pn) ^1 "^X, m^) 

4 



4 \ 

--^2 {q,mx,mn) J 
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7. Diagram D12b 

Here, q = 

a2{m\ - ml) -1 



5 5 
-:^Pk ■ Bi (g, niK, m^) - px ■ Bi {q, m^, uik) - -Pn ■ (g, friK, m^) - p^ ■ Bi {q, m^, tuk) 

-j^Pk ■ Bq (g, rriK, m^) [-2?,m\ - 13m^ + 40A;„ ■ Pk - l^K ■ Ptt + lOpx ■ P-k) 
5 

■ Bq (g, m^-, m^) {-ml + K ■ Pk + K ■ Pw - Pk ■ P-n) 

■ Bq (g, m^, mi^) (-m^ + K ■ Pk + K ■ Pn - Pk ■ Pn) 
-^Ptt ■ 5o (g, rrin, mx) {-I7m\ - ml + 30K ■ Pk) 

--^Pk ■ Bq (g, m^, m/^) (-17m^ + 5m^ + 24A;^ ■ px - QK ■ Ptt + QPk ■ Pn) 

-^P^ ■ Bq (g, m/f, m^) ■ p/^ - 23 (m|. + m^)) 

3 3 
+-k^ ■ Bq (g, mx, m^) ■ /c^ - 3k^ ■ Bq (g, m^, m^) ■ + -k^ ■ Bq (g, m^, m^) ■ fc^r 

3 3 
-3k^ ■ Bq (g, m^, mx) ■ Pn - -PK ■ Bq (g, mK, m^) ■ px - -pK ■ Bq (g, m^, mi^-) ■ px 

-3p^ ■ Bq (g, mx, m^) ■ px - 3p„ ■ Bq (g, m^, mx) ■ Pk 

(-3m^ + ml + 6k„ ■ px - 4:kj, ■ p^, + Apx ■ Ptt) {lOk„ ■pK-'i {m\ + ml)) Bq {q,mx,m^) 

"Tii ^"^"^^ ^"^^ " ^^"^ ' ^^'^ ' ^'^ ~ ' ^'"'^ '^"^^ "^^'^ " ^'^^'^ ' -^^^ ^° 

(l'^"^/^ - 11"^^ - IS^T ■ Pk + 12/c,r ■ Ptt - 12p/< ■ p,r) (l7m^ + - SOA;^ ■ pi^) 5o (g, m^, m^) 

o04 

+- {3mx + ml- 8k^ ■ px + 2k^ ■ - 2px ■ P-k) Bi (g, mx, m„) 
+^ i^'^K + - SA;^ ■ Pa- + 2k^ ■ p^ - 2px ■ Pn) Bi (g, mx, m^) 

io 

-^^2 (g, mx, m„) - \b2 (g, m^, m^-) ) 
18 D / 



20 



8. Diagram D12c 



Here, q=\{pK- k^] 



/c^ ■ Bq {q, rriK, mx) ■ K - '^K ■ Bq (g, uik, mx) ■ Pk + ■ Bq (g, m^, m^) ■ /c^ 
-Ak^ ■ Bo {q, m^, m^) ■ pk - 2p^ ■ Bq {q, uik, m^) ■ px - ■ Bq (g, niK, mx) ■ p^ 
-4p^ ■ Bq (g, m^, m^) ■ - ■ 5o (g, m^, m^) ■ 

+^ (-2m|- - 5ml + 5A;^ ■ p^) {2m% + 3ml + 2k^ ■ px - 3k^ ■ Pn - 2pK ■ Pn) Bo (g, mx, 
+ 7^ {5k^ ■ Pn - 7rnl) [5ml + 2k^ ■ px - 3k^ ■ p^ - 2pK ■ Pn) Bo (g, m^, m^) 

io 

i^K ■ Pn - 4m^) (5ml + 2k„ ■ px - Sk^ ■ Pn - 2pK ■ Pn) Bo (g, m^, m^) 

io 

~T^"^^ ("^"^^ + S^TT ■ Pk - ■ Pn - GpK ■ Pn) Bo (g, m^, m^) 

io 

+- (2m^ + Ami + K ■ Pk - ■ Pn - Pk ■ Pn) Bi (g, mx, mx) 
2 

+- [Qml + k^-px - AK -Pn-PK ■ Pn) Bi (g, m^, m^) 

+^ {9ml + ■ px - 8k^ ■ p„ - 2px ■ Pn) Bi (g, m^, m^) 

1 1 4 \ 

+ (g, m^, m^) - -B2 (g, m^, m^) - -B2 (g, m^, m^) j 



9. Diagram D3 



,7r 



-ik^)\0^'^'^\K\px)n~{p^)) 



a2{m\ - ml) 



-1 



D3 



4/57r2 



k 



-IS/c^ ■ k^ + 13A;^ ■ + IS/c^ ■ + llpx ■ Px + 47pii- ■ p^ + 12p^ • p^) Aq {mx) 



32 

45 
16 

+— (-8A;^ • A;^ + 13A;^ ■ px + ■ Pn + 2pi^ ■ px + 17^^ ■ Pn + ■ p^) Ao {m^) 
y 

16 

+ — {-Ak^ ■ k^ - k^ ■ Px + Ak^ ■ Pn + GPK ■ Px + ^^PK ■ Pn + Pn ■ Pn) ^0 ("^r,) 

15 

544 , , , 112 ^ , , 16 ^ , , \ 
+— {mx) + — (m,) + — K) j 
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10. Diagram Fl 



F3 ji^2 



-y {^aiK ■ Pk - 4a2m|.) Aq {ttik) - ^ [2aik^ ■ px - oi2m\) Aq (m^) 

32 64 \ 

— ^ {^aiK ■ Pk - a-zm]^) Aq (m^) - QAaiAi (niK) + y ai^i (m^) j 



11. Diagram F2 



Here, q = \{K - Pn) 



' F2 ji^2 

6ai/c^ ■ Bi (g, ttik-, m^) - 2aiA;^ ■ Bi (g, m^, m^) + Qaipx ■ -Bi (g, "t,/^, m^) 

3 1 
--tti/c^ ■ Bo (g, m^, m^) + -aik^ ■ Bq (g, m^, m^) - 2aiPA' ■ (g, "^i^, 

1 3 
+ 2«iPa' ■ 5o (g, m^) ml - -aipx ■ Bq (g, rriK, m^) ml 

+ ^aiml (8m|. + 13m^ - 20K ■ Pk + ^Ok„ ■ p^ - lOpx ■ Pn) Bq (g, m^, m^) 

-^aiml (-8m|. + 5m^ + 12k^ ■ px - 6k^ ■ p^ + Gpx ■ Pn) Bq (g, m^, m^) 

+ (-4m|. - 9m^ + ■ p/< - 5k^ ■ p^ + Sp^ ■ Pn) (g, "i^, m^) 

+ {-Am\ + ml + 6kn ■ Pk - ^k^ ■ Pn + ^Pk ■ Pn) Bi (g, mx, m^) 

+ y Qi^s (g, mx, m^) + ^cti^s (g, m^, m^) ) 
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12. Diagram A3 



A3 p^2 

128^2 (-2m^ + ■ Pk + ■ Pn + -ipK ■ Pn) [ml; - ml) Ap (mj^) 

3 {2ml - 2/^7r ■ Pk - 2A;^ ■ p^r + 2pK ■ Pn) 
128a2 (-2m2 + 2k^ ■ pk + 2k^ ■ p^ + ipK ■ Pn) {rn\ - ml) Aq (mx) 



64a2 



3 {2ml - '^^TT ■ Pk - ■ p^ + 2pK ■ p-n) 
-2ml + 2/c7r ■ Pk + 2^ ■ p^ + 4p/< ■ Pn) {mj^ - ml) Aq (m^ 



64^2 



3 {2ml - "^kn ■ Pk - '^K ■ Pn + 2pK ■ Pn) 
-2ml + '^kn ■ Pk + 2/c^ ■ Pn + 4pi^ ■ Pn) {m\ - ml) Aq (m^ 



64ai 



9 {2ml - "^kn ■ Pk - 2/c^ ■ Pn + 2pK ■ Pn) 
-2ml + 2A;^ ■ Pk + 2A;^ ■ Pn + 4p/^ ■ pn) Ai {mx) 



32ai 



3 {2ml - "^kn ■ Pk - '^K ■ Pn + 2pK ■ Pn) 
-2ml + 2A;,r ■ Pk + '^K ■ Pn + 4p/c ■ Pn) Ai (m^) 



32ai 



2ml - '^^■n ■ Pk - '^K ■ Pn + '^Pk ■ Pn 
-2ml + '^^^ ■ Pk + '^K ■ Pn + "^Pk ■ Pn) M {mn) 



3 {2ml - "^^TT ■ Pk - '^K ■ Pn + 2pK ■ Pn) 



APPENDIX B: FINITE VOLUME EFFECTS 



The following argument follows Ref. 12j, where more detailed discussion can be found. 
The relevant correlation function is 



C{t)= [ d^e^-"" {7r{kn)\Hw{0)OKn{t,x)\0) 
Jv 

^VY, {7^{kn)\Hw{0)\K7r, n, P) e~^"* {Ktt, n, P|O^^(0)|0) 



(Bl) 



where contributions from excited state e.g. four particle states are ignored since we are 
interested in asymptotic limit in which t — oo. 

If the volume is sufficiently large, the summation can be approximated with integration : 

C{t) = (B2) 

POO 

V / dE pv{E)e-^' {7r{K)\Hw{0)\K7r,{E,P)) {K7r,{E,P)\OKn{0)\0) , 
Jo 



where the interpretation of pv{E) is given in Ref. [12] • Roughly, py can be understood as 
density of states. 
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Meanwhile, the large volume allows us to rewrite the correlation function, C{t), in terms 
of infinite volume asymptotic states: 



C{t) = (B3) 
(^) / ^^(^(k.)I^H^(O)|i^(pi)vr(p2))(i^(pi)7r(p2)|0...|O). 
We perform a small manipulation on the above equation using 

= {2ti)H{V -pk- p.) {K{pk)ti{p^)\OkM\^) e"""* 

where 



E = Ek{Pk) + K(P.), ^,,(p) = \Jmf + p2. (B4) 
Then the correlation function can be written in terms of 4-momentum integral: 

(7r(k^)|iJH/(0)OA'.(t,P)|0) = 

{n{K)\HwmK{pKMp^)){K{pK)Ap.)\OKM\^) 

where 

Pi = {E,,pi). (B5) 

Since the integral is in covariant form, one can easily change the integration variables with 
Lorentz transformed ones which brings the two particle state into CM frame: 

(7r(k^)|ffw'(0)Ox.(t,P)|0) = 

^2^)' / / (^^''(^^ - O^^i^l - <nE* -E],- E:)5\c,k + q.)e-^* 

Note that the amplitudes are treated as scalars. Using the (5-function, one can do some 
integrations, 

{7r{K)\Hw{0)OKAt,Pm = 

(2.)- / dE e- / ^^^(^^^(i^* - EM - EM)) 

(7r(k;)|ifH'(0)|K(q)7r(-q)) (K(q)7r(-q)|Ox.(0)|0) , 
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where E]^ and E* are fixed respectively by 

EM = 



E^(q) = v/m2+q2. (B6) 

In order to proceed, we have to assume that the creation operator Ok-k is chosen so 
that (J'r(q)7r(— q)|(9x^(0)|0) does not have angular dependence. In fact, the contribution 
originating from the particular form of Oxn will be canceled in the end. The role of Oxn is 
restricted to keeping the initial states in S-wave. With this assumption, a further simplifi- 
cation is possible : 



where q* is defined as 



dE e 



-Et 



1 E 



q* EK{q*) + E^{q 
dn {7r{K)\HwmK{cL)7r{-cO) 



(J^(q*)7r(-q*)|O^.(0)|0) 



(B7) 



E* 



nil 



-r g'" + V + ?*^- (B8) 
Since (ii'(q*)7r(— q*)|Ci^^(0)|0) does not have any angular dependence, we write this as 
(/s:7r,E*|Ox7r(0)|0). Moreover, / (7r(k;)|i7vF(0)|i^(q)7r(-q)) is also a function of E* 
because of angular averaging so we use a notation dVlE*- 
By comparing Eq. flB2|) and Eq. flB7l) . one can deduce 

Pv{E) (7r(kJ|i7H^(0)|K7r, (E,P)) (AV, {E,V)\OkM = (B9) 
(J^) (i^7r,E*|O^.(0)|0) j dQE' {7r{k:)\H^mK{q)7r{-cD) . 

However, Eq.( ]B9l) is not satisfactory because of the appearance of an unknown factor 
(A'tt, i?*|(9j^7r(0)|0). In order to eliminate this factor, we consider a correlation function, 

'OKAt,P)OL{0,P)) . (BIO) 

One can easily imagine that by the similar argument, the following equation can be obtained. 



TT 



(2vr) 



^^]\{Kn,E*\OKA0)\0)\'. 



(Bll) 



By dividing Eq. flBQP by the square root of Eq. flBllI) . the relation between matrix elements 
can be deduced, 



{7r{K)\HwmKrt,{E,P)) 
1 1 



TT 



47r Vp^V (27r)3 V^* 



dnE* (7r(k;)|i7H/(0)|/^(q)7r(-q)) 
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APPENDIX C: LEGS FOR Of^J^J/^ AND O^^^'^'"^ 



It has been argued that LECs for operators with (8,8) AJ = 1/2 can determined without 
K 111^1=0 in Ref ■ [lo| . In this section, we will show that LECs for operators with (27,1) 
A/ = 1/2 as well can be determined without calculating A/ = 1/2, K — > 71171=0 matrix 
elements. Its analytic term at physical kinematics is given by 
4:i{mj^ — m^) 



P 



+2 {-6di - 24 + 24 + 64 + - 2^20 + 8^24) ml 

The K ^ TT transition amplitude with (27, 1) A/ = 3/2, are : 
1 



(CI) 



4:a27{PK ■ K) + 64^24 (Pi^ ■ /^tt) +8 {d^ + dj - ^20) "^^(pi^ ■ K) 



+m% (8 (^4 + 2^7 - ^20) (Pe- ■ K) - lQd2ml) 



(C2) 



From the above, one can determine 027, d2, d^, ^20 — (^4, (^24 by varying masses and momenta. 

The xPT formula for K 7^7^1=2 transition amplitude with unphysical kinematics cor- 
responding to the initial kaon and the final pions are at rest, is given by 



ot27{fnK + 'mT,)/2 + (^4 + rfs + Ad-j - d2o) m 



K 



P 

f 2 (ci2o - ^2) nxT^m]^ + (8^4 + ci5 + 2^7 - 3c?2o) mlniK + 2^2^^ 



(C3) 



By inspecting the above formula, we can see that d^ — d^ can be determined. Since the xPT 
formula for the K — > vac transition amplitude is 

4:8idi {m\ — m^)^ 



/ 



(C4) 



di can be fixed from this. Finally, d^ can be determined by using A/ = 1/2 A' ^ tt transition 
amplitude whose xPT formula is 
1 



P 



4a27(px ■ K) + (8 {di - Sde + 2^7 - ^20) (pk ■ K) - 16 (3c/i + ^2) "2^) 
+48c?im^ + 64:d2i{pK ■ + 8 (^4 + Mq + dj - d2o) {px ■ K)^^ . 

(C5) 

Since ^27, di, ^2, d^—d^, dg, dj, ^20— c?4, (^24 are determined, one can reconstruct the physical 
K ^ 7171 transition amplitude with (27, 1) A/ = 1/2 using Eq. OCip . 
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